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Dynamic Analysis of Geared Rotors 
by Finite Elements 


A finite element model of a geared rotor system on flexible bearings has been 

shafts' n he ?n eI m . Cl “ deS the rotar y inertia of shaft elements, the Lai loading 
onshafls, flexibility and damping of bearings, material damping of shafts and the 

ofgear mes \ The coupling between ,he lorsional and 
transverse vibrations of gears were considered in the model. A constant mesh stiff- 
ness was assumed The analysis procedure can be used for forced vibration analysis 
f geared rotors by calculating the critical speeds and determining the response of 

JiLf ° n , , hS Sha/ ' S - ,0 m0SS unbalances > geometric eccentricities of gears, and 
displacement transmission error excitation at the mesh point. The dynamic mesh 
forces due to these excitations can also be calculated. The model has been applied to 
several systems for the demonstration of its accuracy and for studying the effect o f 
bearing compliances on system dynamics. . 8 11 1 


Introduction 

Even though there have been numerous studies on both 
rotor dynamics and gear dynamics, the studies on geared rotor 
dynamics have been rather recent. The study of the dynamic 
behavior of geared rotor systems usually requires the coupling 
between torsional and transverse vibration modes to be includ- 
ed m the model, although this is not a problem for rotors 
without gears. 

Although several modeling and solution techniques such as 
lumped mass models and the use of the transfer matrix 
method have been applied to rotor dynamics problems, the 
finite element method seems to be a highly efficient way for 
modeling such systems. In one of the early examples of the 
finite element method, Nelson and McVaugh (1976) used a 
Rayleigh beam finite element including the effects of the 
translational and rotary inertias, gyroscopic moments, and the 
axial load. The work of Zorzi and Nelson (1977) was the 
generalization of the previous study (Nelson and McVaugh, 
1976) to include internal damping. Later, Nelson (1980) 
developed a Timoshenko beam by adding shear deformation 
to the Rayleigh beam theory. This model was extended by 
Ozguven and Ozkan (1983) to include effects such as 
transverse and rotary inertia, gyroscopic moments, axial load, 
internal hysteretic and viscous damping, and shear deforma- 
tions in a single model. None of the models described above 
can handle geared rotor systems, although they are capable of 
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determining the dynamic behavior of rotors which consist of a 
shaft supported at several points and carrying rigid disks at 
several locations. 


Oear dynamics studies, on the other hand, have usually 
neglected the lateral vibrations of the shafts and bearings, and 
have typically represented the system with a torsional model 
Although neglecting lateral vibrations might be a good ap- 
proximation for systems having shafts with small com- 
phances, it was observed experimentally (Mitchell and Mellen, 
, .) 1 , dynamic coupling between the transverse and 

torsional vibrations due to the gear mesh affects the system 
behavior considerably when the shafts have high compliances. 
This fact directed the attention of investigators to the inclu- 
sion of lateral vibrations of the shafts and the bearings in 
mathematical models. Lund (1978) developed influence coeffi- 
cients at each gear mesh by using the Holzer method for tor- 
sional vibrations and the Myklestad-Prohl method for lateral 
vibrations. He obtained critical speeds and a forced vibration 
response by coupling the results of these two methods. Hamad 
and Seireg (1980) studied the whirling of geared rotor systems 
supported on hydrodynamic bearings. Torsional vibrations 
were not considered in this model and the driven shaft was 
assumed to be rigid. Iida et al. (1980) considered the same 
problem by taking one of the shafts to be rigid and neglecting 
the compliance of the gear mesh, and obtained a three degree 
of freedom model to determine the first three vibration modes 
and the forced vibration response due to the unbalance and 
the geometric eccentricity of one of the gears. They also 
reported that their theoretical results confirmed experimental 
measurements. Later, Iida et al. (1984, 1985, 1986) applied 
their model to a larger system which consists of three shafts 
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coupled by two gear meshes. Hagiwara et al. (1981) developed 
a simple model that includes the transverse flexibilities of the 
shafts by using discrete stiffness values, and studied the forced 
response of geared shafts due to unbalances and runout er- 
rors. They took the damping and compliances of the journal 
bearings into account and assumed the mesh stiffness to be 
constant. 

Another approach in the study of gear dynamics has been to 
use transfer matrix methods. In most of these studies the mesh 
stiffness was taken to be periodic. Daws (1979) developed a 
three-dimensional model considering the mesh stiffness as a 
time varying three-dimensional stiffness tensor. He included 
the force coupling due the interaction of gear deflection and 
time varying stiffness, whereas, he neglected the dynamic 
coupling in the model. As a continuation of this study, 
Mitchell and David (1985) showed that dynamic coupling 
terms were dominant on the dynamic behavior of the system. 
Iwatsubo et al. (1984a) also used the transfer matrix method in 
their models and calculated the force response due to mass un- 
balance only for a constant mesh stiffness. Later, they (1984b) 
included the effects of periodic variation of mesh stiffness and 
profile errors of both gears. 

Neriya et al. (1984) extended the model of Iida et al. (1980) 
by representing a single gear by a two mass, two spring, two 
damper system which used a constant mesh stiffness. The 
shafts were assumed massless and equivalent values for the 
lateral and torsional stiffnesses of shafts were used to obtain a 
discrete model. Later Neriya et al. (1985) employed the finite 
element method in finding the dynamic behavior of geared 
rotors. They also found the forced vibration response of the 
system due to mass unbalances and runout errors of the gears 
by using modal summation. Bagci and Rajavenkateswaran 
(1987) used spatial finite line-element technique to perform 
mode shape and frequency analysis of coupled torsional, flex- 
ural and longitudinal vibratory systems with special applica- 
tion to multicylinder engines. They concluded that coupled 
torsional and flexural modal analysis is the best procedure to 
find natural frequencies and corresponding mode shapes. An 
extensive survey of mathematical models used in gear 
dynamics analyses is given in a recent paper by Ozguven and 
Houser (1988a). 



Fig. 1 A typical gear-rotor system 


The objective of this study was to develop a finite element 
program for the dynamic analysis of geared rotor systems and 
to study the effect of bearing flexibility, which is usually 
neglected in simpler gear dynamics models, on the dynamics 
of the system. In the formulation of rotor elements, except for 
gears, the rotor dynamics program ROT-VIB which has been 
developed in a previous study (Ozguven and Ozkan, 1983; 
Ozkan, 1983) was used. However, due to the coupling between 
torsional and transverse vibration modes, a torsional degree of 
freedom has been added to the formulation, and some special 
features of ROT-VIB have been omitted. 

Theory 

A typical geared rotor system is shown in Fig. 1. The system 
consists of a motor connected to one of the shafts by a cou- 
pling, a load at the other end of the other shaft, and a gear 
pair which couples the shafts. Both shafts are supported at 
several locations by bearings. Hence, a geared rotor system 
consists of the following elements: 

(a) shafts, 

( b ) rigid disks, 

(c) flexible bearings, _ 

( d ) gears. 

When two shafts are not coupled, each gear can be modeled as 



Nomenclature 


[C] = damping matrix of the 
system 

. c m = mesh damping coefficient 
[C m j = mesh damping matrix 
c s = modal damping value of 
5th mode ‘ • 

d ]t d 2 = diameters of driving and 
driven shafts,, respectively 
e g , e p = geometric eccentricities 
of the gear and the pin- 
ion, respectively 
e t = amplitude of the static 
transmission error 
excitation 

F s = average value of force 
transmitted (static load) 
[F, ) = total force vector of the 
system 

/ = unit imaginary number 
I g y l p = mass moment of inertias 
of the gear and the pin- 
ion, respectively 

J d , J m = mass moment of inertias 
of load and motor, 
respectively 


k m = mesh stiffness coefficient 
[K m ] = mesh stiffness matrix 
K td = torsional compliance of 
the flexible coupling 
kxx> kyy - bearing stiffness values 
Lj, l 2 = lengths of driving and 

driven shafts, respectively 
m g , m p = masses of the gear and 
the pinion, respectively 
A/,, M 2 = moments due to the 
dynamic mesh forces 
N p = tooth number of driving 
gear 

[ q ) = total response of the 
system 

r gy r p = base circle radii of the 
gear and the pinion, 
respectively 
t = time 

U gi Up = mass unbalances of the 
gear and the pinion, 
respectively 

W x , W 2 = dynamic mesh forces 
x gt x p = coordinates perpendicular 


to the pressure line at the 
centers of the gear and 
the pinion, respectively 
y gt y p = coordinates in the direc- 
tion of the pressure line 
at the centers of the gear 
and the pinion, 
respectively 

[a] ~ dynamic compliance 
matrix 

0 2 , = total angular rotations of 

the gear and the pinion, 
respectively 

6 p , 0 g = fluctuating parts of 0, 
and^0 2 , respectively 
o)p , o) g = rotational speeds of driv- 
ing and driven shafts, 
respectively 

o) r = rth natural frequency 
[$] = modal matrix 
[4> J ) = 5th normalized eigen 
vector 

= modal damping coeffi- 
cient for the 5th mode 
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SwrtS as’ ‘ h ' dynam “ m " h f “'“ ln Erection 

c m (y p + f p 6 1 + e p u p cosd ] —y g —r g 0 2 — e g u g cosd 2 
-e.N^pCo^Npd^)) 

+ k m(y p + r p 6 x +^11101 ~y g — r g d 2 -e g sind 2 -e,sin(N p d t )); 


iv 2 =-yv xl 


a rigid disk. However, when they are in mesh, these rigid disks 

- C0 "r ed by a sprin g- damper element representing the 
mesh stiffness and damping. 

For the formulation of the first three types of elements 

oSn i°o V 8 e ,\ * eXiSt l ng pr0gram ROT - VIB . (Ozguven and 
Ozkan, 1983) was used. ROT-VIB is a general purpose rotor 

dynamics program which can calculate whirl speeds, cor- 
respondmg mode shapes, and the unbalance response of shaft- 
rigid disk-bearing systems by including the effects of rotary 
and transverse inertia, shear deformations, internal hysteretic 

ratvid US u’ an ) pin 8 , 3X131 !°ad and gyroscopic moments. In 
KOT-VIB, the classical linearized model with eight spring and 
amping coefficients are used for modeling bearings, and 
finite elements with four degrees of freedom at each node (ex- 
cluding axial motion and torsional rotation) are employed for 
the shaft elements. 

In the present analysis, the formulation used in ROT-VIB 
for these elements was employed with some modifications, 
rirst, m order to avoid nonsymmetric system matrices which 
result in a complex eigenvalue problem, the gyroscopic mo- 
ment effect was ignored and the internal damping of the shafts 
was included only in the damping matrix. Second, the gear 
mesh causes coupling between the torsional and transverse 
vibrations of the system, which makes it necessary to include 
the torsional degree of freedom. Therefore, the mass and stiff- 
ness matrices of the system which are taken from ROT-VIB 
have been expanded in the new model to include the torsional 
motion of the shafts. Hence, five degrees of freedom have 
been defined at each node with only axial motion being ex- 
cluded. The axial motion, which would be important for 
helical gears, could easily be included in subsequent analyses. 

Gear Mesh Formulation. A typical gear mesh can be 
represented by a pair of rigid disks connected by a spring and a 
damper along the pressure line which is tangent to the base 
circles of the gears as shown in Fig. 2. In this model, the mesh 
stiffness and damping values are assumed to be constant. The 
tooth separation is not considered since the gears are assumed 
to be heavily loaded. By choosing the y axis on the pressure 
line and the at axis perpendicular to the pressure line, the 
transverse vibrations in the x direction are uncoupled from 
both the torsional vibrations and the transverse vibrations in 
the y direction, provided negligible mesh friction. For the 
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0) 

, ( 2 ) 
where W , and W 2 are mesh forces in th ty n and y directions 
at the pinion and the gear locations respectively. Here, c m and 
m are mesh damping and mesh stiffness values, e„ and e are 
f-ZT ecce " tric ' ues of the p mion and the gear, and r* and 
and fl 5 C th C 6 ^ 1 “ ° f thC Pmi0n 3nd the gear - The an § les 8, 

and 0 2 are he total angular rotations of the pinion and the 
gear, respectively, and are equal to 

°i = 9 p +a p t; e 2 =d g + u gt (3,4) 

where 6 p and 10, are the alternating parts of rotations and w 
nd o) g are the spin speeds of the driving and driven shafts^ 
respectively. The displacement ,, which may be considered to 

nnin, nTfr err ° r ? xcitation - is a PPKed at the mesh 
point. This displacement is usually taken to be sinusoidal at 

the gear mesh frequency, but one could include higher har- 

fl988bl thm i/ 1 ^ b vfi en Sh ° Wn by ° 2 8 uven and Houser 
(1988b) that it is possible to simulate the variable mesh stiff- 
ness approximately, by using a constant mesh stiffness with a 
displacement excitation representing the “loaded” static 

t^rir:- Thus ’ by choosin 8 e, as the amplitude of 
the loaded static transmission error, the effect of the variable 

. . h S “ f / ness can be approximately considered in the model. 
Mesh forces also cause moments about dynamic centers of 
the gears which are equal to 

M l = W l (r p + e p cosd l ); M 2 =W 2 (r g +e g cos0 2 ) (5,6) 

^ th Vr al angalar P osi tion of geometric eccentricities 
l . ke " t0 be zero - Tbe mesh stiffness and damping matrices 
and the force vector of the system due to gear errors and un- 
balances can be obtained by writing the force transmitted as 
the summation of the average transmitted force (static load), 
F s and a fluctuating component, and then neglecting high 
order terms following the substitution of equations (1) and (2) 
into equations (5) and (6). Defining the degrees of freedom of 
the system at which the gear coupling effect appears, as 

(Qi } = 0 ^, o p y g e g ] T (7) 

the additional mesh stiffness matrix which causes the coupling 
effect and corre spon ds to (q, ) can be obtained from equa- 
tions (1), (2), (5), and (6) to be 

k m 




k n,r p 

-k„ 


k m r p 

k m>% 

~ k n,r p 


— k„ 


~ k n,r g 


~ k m r p k m r p 


k m r t k m r p r g k m r g 


k n,r g 
k ^g J 


( 8 ) 


Similarly, the mesh damping matrix can be found to be 


[CJ = 


c 

*~m r p 




C m r p C m r p r g 




C m r g C m r p r g 


C m' 2 g J 


(9) 


The other degrees of freedom defined at nodes p and g have 
not been included in the vector {q, ) since elements of [K m ] 
and corresponding to these degrees of freedom are all 
zero. For the degrees of freedom expressed as 
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Table 1 Parameters of the system shown In Fig. 1 


J m : 0.459 Kg-m 2 
J d : 0.549 Kg-m 2 
Ip : 0.030 Kg-m 2 
Ig : 6.28x10 3 Kg-m 2 
m p : 16.96 Kg 
m g : 5.65 Kg 
r p : 0.1015 m 
r g : 0.0564 m 
L, : 0.78 m 
L 2 : 0.40 m 


dj : 0.03 m 

d 2 : 0.02 m 

e g 1.2xl0’ 5 m 

U g : 2.8xl0’ 4 Kg-m 

e p : 0.0 

U p : 0.0 

e t : 0.0 

K (C j : 115.0 N-m/rad 

k m : 2.0x1 0 8 N/m 


lq 2 ) = b >p x p e p y g x g $ g ] T (io) 

the force vector due to the gear runouts, the static transmis- 
sion error, and the gear mass unbalances U D and U g is given 
by 

U p o)ls\nu p t + Ft 
U p w 2 p cosw p t 


[FJ = 


-F s e p cosw p t + r p F x 
UgO) 2 g s\nu g t-F l 


U g UgCOS(j) g t 

^ F s e g cosw g t-rgF x y 


( 11 ) 


where 

F\ =c m (egWgCOso>gt-e p o) p cosw p t + e t N p G) p cos(N p o) p t)) 

+ k m (e g sino) g t - e p sin w p t + e t s\n(N p o) p t) ) (12) 

Adding the mesh stiffness matrix given by equation (8) to 
the stiffness matrix of the uncoupled rotor system yields the 
total stiffness matrix of the system. The natural frequencies u r 
and the mode shapes (u r ) of the system can be determined by 
solving the governing eigenvalue problem. In the solution, the 
Sequential Threshold Jacobi method was used. 


Forced Response. The total force vector is obtained by 
combining the force vector due to the mass unbalances of the 
shafts and the other disks and the force vector due to the mass 
unbalances of gears and gear errors as given in equation (11). 
This vector is the sum of harmonic components with three dif- 
ferent frequencies a> p , cj g and (N p u p ) , and has the following 
general form 

t F J = IF^Jsinc o p t+ {F cp ) cosay + {F sg )sinay+ [F C J cosay 

+ { F, m ) sin (N p u p t) + { F cm ) cos (N p u p t) (13) 

The total response of the system due to this excitation can be 
written as 


f Q 1 = [oc p ] f F sp ) sinoy + [a p ] [ F cp ) cosay 

+ [a g ] l F sg Jsina y + [a g ] ( F cg ) cosay 
+ [«ml f *sm ) sin ( N p <a p t) + [a m ] [ F cm ) cos (N p o> p t) (14) 

where [a p ], [a g ] and [a m ] are the dynamic compliance 
matrices corresponding to the excitation frequencies, oj p , u g 
and (N p u p ), respectively, and given by 


ki= L- 


s= 1 




05) 


ki=D 


I* 1 )!* 1 ) 


" (0)|-u| + /(0 J C J ) 


(16) 



Rotational Speed of Shaft 2 (Hz) 

Fig. 3 Comparison of the theoretical values of the dynamic deflection 
of the driven gear in the pressure line direction with the experimental 
results given by lida et al. (1980) 



WU* 5 ) 7 

(wl-Nlul + iNpWpC;) 


(17) 


Here, { & ) represents the 5th mass matrix normalized modal 
vector, n is the total number of the degrees of freedom of the 
system, / is the unit imaginary number, and c s is the 5 th modal 
damping value given by the^yth diagonal element of the 
transformed damping matrix [C] where 

lC\ = m T lQm (18) 

and [4>] is the normalized modal matrix. In this approach it is 
assumed that the damping matrix is the proportional type, 
which is usually not correct for such systems. When the damp- 
ing is not proportional, the transformed damping matrix [C] 
will not be diagonal in which case c s will still be the 5 th 
diagonal element and all nonzero off-diagonal elements are 
simply ignored when using the classical uncoupled mode 
superposition method. Another approach for including damp- 
ing in the dynamic analysis of such systems would be assuming 
a modal damping, for each mode and then replace c s in 
equations (15) to (17) by 2f 5 oj s . However, it is believed that us- 
ing the actual values for damping, when they are known, and 
employing an approximate solution technique may result in 
more realistic predictions than assuming a modal damping 
value for each mode. 


Applications and Numerical Results 

Comparison with an Experimental Study. As the first ap- 
plication, the experimental set-up of lida et al. (1980) has been 
modeled. As shown in Fig. 1, the system consists of two 
geared rotors, one is connected to a motor with a mass mo- 
ment of inertia of J m and the other is connected to a load with 
a mass moment of inertia of J d . Each shaft is supported by a 
pair of ball bearings. The parameters of the system are listed 
in Table 1 . The gears with inertias I p and I g are both mounted 
on the middle of the shafts of lengths L { and L 2 and diameters 
d x and d 2i respectively. In their study, lida et al. (1980) did not 
specify the length of the second shaft, L 2 , and the properties 
of bearings and couplings. Instead they gave the total tor- 
sional stiffness values for driving and driven parts of the 
system and a total transverse stiffness value for the second 
shaft. Therefore, the length of the second shaft, 1 2 , and the 
torsional stiffness of the first coupling, K tc 1, have been 
estimated so that the total values given by lida et al. (1980) 
were obtained. The forced vibration response due to the 
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Table 2 Parameters of the system shown in Fig. 4 


J m : U5xl O' 2 Kg-m 2 
J d : 5.75x1 0‘ 3 Kg-m 2 
m m * 9.2 Kg 
m d : 4.6 Kg 
I ? : 1.15xIO" 3 Kg-m 2 
I 8 : 1.15xl0' 3 Kg-m 2 
F s : 2500 N 


r p : 0.047 m 
r g : 0.047 m 
m p : 0.92 Kg 
m g : 0.92 Kg 
k xx , k yy : variable 
k m : 2.0x1 0 8 N/m 


Table 3 Parameters of the system shown in Fig. 6 


Ip : 0.0018 Kg-m 2 

ITIp 

1.84 Kg 

I g : 0.0018 Kg-m 2 

m g 

1.84 Kg 

r p : 0.0445 m 

c t 

9.3xl0’ 6 m 

r g : 0.0445 m 


1. Ox 10 s N/m 

kxx» kyy : variable 

N p : 

28 teeth 


geometric eccentricity e g and the mass unbalance U g (e n e py 
and U p are all zero) is shown in Fig. 3, along with the ex- 
perimental results of Iida et al. (1980). Since no information is 
given about the damping values of the system, a modal damp- 
ing of 0.02 has been used at each mode in the computations. 
As seen in Fig. 3, predictions from the analytical model show 
good correlation with the experimental results. 

Response due to Geometric Eccentricities, Mass Un- 
balances, Static Transmission Error, and Mesh Stiffness 
Variation. As a second example, the system which was used 
by Neriya et al. (1985) was studied to investigate the effects of 
geometric eccentricities and mass unbalances of the gears on 
the forced response of the system. The natural frequencies, 
mode shapes, and the responses at both gear locations due to 
geometric eccentricities and mass unbalances of gears were 
almost identical to those which were documented by Neriya et 
al. (1985). The results of this analysis have not been included 
in this study since the gear eccentricities and unbalances excite 
the system at the shaft rotational frequencies as was shown in 
the first example. The contribution of such low frequency ex- 
citations on the generated gear noise is usually negligible when 
compared with those of high frequency excitations caused by 
the static transmission error and the mesh stiffness variations. 

On the other hand, the system shown in Fig. 4 has been 
modeled to obtain the dynamic mesh force due to the static 
transmission error excitation of amplitude e t and frequency 
N p o) p representing the mesh stiffness variation. Dimensions of 
the rotors are shown in Fig. 4 and other system parameters are 
listed in Table 2. The bearings are assumed to be identical and 
geometric eccentricities and mass unbalances for gears are 
assumed to be zero, so that only excitation causing a forced 
response is the static transmission error excitation defined. 
Since the displacement input approximates the loaded static 
transmission error, the value of e , was taken as the amplitude 
of the loaded static transmission error. Figure 5 shows the 
variation of the ratio of dynamic to static mesh load with rota- 
tional speed for three different bearing compliances. The first 
two small peaks of Fig. 5 correspond to torsional modes of 
shafts and the third peak corresponds to the coupled 
lateral/torsional mode governed by the gear mesh. 

As shown in Fig. 5, when the bearing stiffnesses are de- 
creased, the dynamic force also decreases considerably 
because of a resulting decrease in the relative angular rotations 
of the two gears. Although the displacements in the y direction 
increase slightly, they do not appreciably affect the dynamic 
force. In this example, a mesh damping corresponding to a 



Fig. 4 The system of the second example (dimensions are in 
millimeters) 



Rotational Speed (Hz) 

Fig. 5 Variation of dynamic to static load ratio with frequency for three 
different bearing compliances 
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Fig. 6 The system analyzed as a third example (dimensions are in 
millimeters) 


modal damping of 0.1 in the mode of gear mesh has been 
used. This was the value employed by several investigators for 
the same problem. 

The Effect of Bearing Compliances on Gear Dynamics. A 
parametric study on the system shown in Fig. 6 was per- 
formed. The effects of bearing compliances on the natural fre- 
quencies and the forced response of the system to the har- 
monic excitation respresenting the static transmission error 
and the mesh stiffness variation were studied. The system 
parameters are given in Table 3. The natural frequencies and 
the physical descriptions of the corresponding modes for a 
value of bearing stiffnesses k xx = k yy = 1.0 xlO 9 N/m are 
presented in Table 4. The forced response at the pinion Ioca- 
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Table 4 First 14 natural frequencies of the system of Fig. 6 for the case 
■of = 10 . 


Natural Frequency 
(Hz.) 

Description of Corresponding 
Mode 

0 

torsional rigid body 

581 

transverse, torsional 

687 

transverse, x dir., driving shaft 

689 

transverse, y dir. 

691 

transverse, x dir., driven shaft 

2524 

transverse, torsional 

3387 

transverse, y dir. 

3387 

transverse, x dir., driving shaft 

3421 

transverse, x dir., driven shaft 

3421 

transverse, y dir. 

6447 

torsional, driving shaft 

6539 

torsional, driven shaft 

6831 

transverse, x dir., driving shaft 

6840 

transverse, y dir. 

w 0.00007 T T — - T 


< 0.00006 - 
o 

c 0.00005 - 
J 

~ 0.00004 - 

© 

| 0.00003 - 

u 

u 

G 0.00002 


0.00001 


iiij 


0.00000 L 



- k„/k m - 0,1 

- k„/k m = 0.5 
-k„/k m « 1.0 
-k„/k w » 10.0 



1500 2000 2500 3000 


Excitation Frequency (Hz) 

Fig. 7 Forced response of the system shown In Fig. 6 to the displace- 
ment excitation at the direction of pressure line (at pinion location) for 
four different bearing compliances 


tion in both the transverse (pressure line) and torsional direc- 
tions, and dynamic mesh forces are plotted in Figs. 7, 8, and 9, 
respectively. As shown in Figs. 7 and 8, the system has peak 
responses only at two natural frequencies within the frequency 
range considered. Mode shapes corresponding to these two 
natural frequencies are presented in Fig. 10. If the free vibra- 
tion characteristics of these two modes are investigated in 
detail, it is seen that the dynamic coupling between the 
transverse and torsional vibrations at these two modes are 
dominant. It is also seen that dynamic loads are high at only 
the second one of these two modes as shown in Fig. 9. The 
reason for this is that the transverse and torsional vibrations 
for the second mode considered are in phase. This results in 
large relative deflections at the mesh point which implies that 
this mode is governed by gear mesh. It is also seen from these 
Figures that lowering the values of bearing stiffnesses causes a 
decrease in both the values of the natural frequencies and the 
amplitudes of the peak responses and dynamic loads. 

Figure 1 1 shows the variation of these natural frequencies 
with bearing stiffnesses for three different shaft compliances: 

(a) long shafts (low stiffness) with dimensions given in Fig. 6, 

( b ) moderately compliant shafts with half the length of case 



500 1 000 1 500 2000 2500 3000 


Excitation Frequency (Hz) 

Fig. 8 Forced response of system shown in Fig. 6 to the displacement 
excitation at the torsional direction (at pinion location) for four different 
bearing compliances 



Excitation Frequency (Hz) 

Fig. 9 Dynamic to static load ratios for system shown in Fig. 6 due to 
the displacement excitation for four different bearing compliances 




Fig. 10 Mode shapes corresponding to natural frequencies at which 
the system shown In Fig. 6 has peak responses 
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(a) 



(b) 



Bearing Stiffness Factor, k xx /k m 

Fig. 11 Variation of the first and the second natural frequencies con- 
sidered with bearing stiffnesses for three different shaft compliances 


a, (c) very short (stiff) shafts. The shaft and the bearings sup- 
porting the gears can be thought of as two springs connected 
in series. When one of these components is very stiff com- 
pared to the other, its effect on the overall dynamic behavior 
becomes negligible. When the mode shapes for these two 
modes are examined for the case of short shafts and stiff bear- 
ings, the first of these two modes becomes purely torsional, 
while the lateral vibrations become more important in the sec- 
ond mode. As shown in Fig. 11(g), since the mode considered 
becomes purely torsional in the case of short shafts and stiff 
bearings, the value of this natural frequency does not change 
as bearing stiffnesses exceed a limiting value. For the other 
mode considered, the natural frequency becomes very high 
when a short shaft is used with a very stiff bearing, since the 
lateral vibrations are more dominant than torsional vibrations 
in this mode as shown in Fig. 11(6). Similarly, when the 
shafts are flexible enough the effect of bearing stiffnesses on 
the natural frequency becomes negligible above a limiting 
value of bearing stiffness. 

Conclusion 

In this study, a finite element model to investigate the 
dynamic behavior of geared rotor systems has been developed. 
In the analysis, transverse and torsional vibrations of the 
shafts and the transverse vibrations of the bearings have been 
considered. The gear mesh was modeled by a pair of rigid 
disks connected by a spring and a damper with constant values 


which represent average mesh values. Tooth separation was 
not considered. The model developed yields the natural fre- 
quencies, corresponding mode shapes, and forced response of 
the system to the mass unbalances and the geometric eccen- 
tricities of gears and the transmission error. 

Although a constant mesh stiffness was assumed, the self- 
excitation effect of a time-varying gear mesh was included into 
the analysis by using a displacement excitation representing 
the static transmission error. It may be justified to solve time- 
varying equations in simpler models with a few degrees of 
freedom. However, for large models such as the ones used in 
this study, avoiding time-varying stiffnesses and transient 
solutions saves considerable computation time. In the example 
problems only the first harmonic of the static transmission er- 
ror was considered. A good correlation between the predic- 
tions and previous experimental and theoretical studies has 
been found. 

Finally, it has been shown that the bearing compliances can 
greatly affect the dynamics of geared systems. Decreasing the 
stiffness values of bearings beyond a certain value lowers the 
natural frequency governed by the gear mesh considerably. 
However, in the case of compliant shafts, when the bearing 
stiffnesses are above a certain value the natural frequency cor- 
responding to the gear mesh does not change considerably by 
increasing bearing stiffnesses. On the other hand, it has been 
seen that the amplitudes of dynamic to static load ratio and 
the deflections at the torsional and transverse directions are 
decreased by using bearings with higher compliances, which 
suggests that the bearing compliance may also affect the 
dynamic tooth load, depending upon the relative compliances 
of the other elements in the system. 
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